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Abstract 

The hypercomputers compute functions or numbers, 
or more generally solve problems or carry out tasks, that 
cannot be computed or solved by a Turing machine. Sev- 
eral numerical simulations of a possible hypercomputa- 
tional algorithm based on quantum computations pre- 
viously constructed by the authors are presented. The 
hypercomputability of our algorithm is based on the 
fact that this algorithm could solve a classically non- 
computable decision problem, the Hilbert's tenth prob- 
lem. The numerical simulations were realized for three 
types of Diophantine equations: with and without solu- 
tions in non-negative integers, and without solutions by 
way of various traditional mathematical packages. 

1 Introduction 

Although there have been several conferences and var- 
ious specialized publications dedicated to the topic of 
hypercomputation, this is still an area in its gestation 
and development process. The term 'hypercomputer' 
denotes any data processing device (theoretical, poten- 
cially realizable or that can be implemented) capable of 
carrying out tasks that cannot be performed by a Tur- 
ing machine (TM) |1 2|. Proposed models of hyper- 
computation surge from diverse disciplines, such as non- 
classical logics, computability, neural networks, or rela- 
tivistic physics Q][3l|4). Notwithstanding the prolifera- 
tion of theoretical hypercomputation models, the possi- 
bility of a conceivable construction of a hypermachine, 
continues to be controversial and held under study and 
analysis. 

This work presents several numerical simulations of a 
possible hypercomputational algorithm based on quan- 
tum computations previously constructed by the authors 
151 . The hypercomputability of our algorithm is based 
on the fact that this algorithm could solve a TM non- 
computable decision problem, Hilbert's tenth problem 
1 6 1 . Our algorithm is based on the algorithm proposed by 
Tien D. Kieu |7]|8][9), with the difference that we have 
selected the Infinite Square Well (ISW) instead of the 
(one-dimensional) Simple Harmonic Oscillator (SHO) 



as the underlying physical system. Since our model 
exploits the quantum adiabatic process for unbounded 
operators |10|, and the characteristics of the infinite- 
dimensional irreducible representation of the dynamical 
algebra su(l, 1) associated to the ISW 1 1 1 1, the simula- 
tion is only able to partially illustrate our algorithm. In 
other words, this work does not pretend to affirm that 
the numerical simulations constitute our hypercomputa- 
tional quantum algorithm. 

2 The Hypercomputational Quantum Algorithm 

The Hamiltonian operator H and the energy levels 
E n , for a particle with mass m in an ISW with length 
ira are II II 



H = i 
E 



h 2 d 2 



2m dx 2 2ma 2 ' 
(h 2 /2ma 2 )n{n + 2), 



where the action of H on its eigenvectors 

{|n)|neN}, N = {0,1,...}, 



(1) 
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is given by H \ n) = E n | n). Due to the spectral struc- 
ture of the ISW, the dynamical algebra associated with 
it, is the Lie algebra su(l, 1) II II . This is a tridimen- 
sional algebra that satisfies the commutation relations 
[K-,K+] = K 3 and [K ± ,K 3 ] = T^K ± , where op- 
erators K + , and K 3 are called creation, annihilation 
and Cartan operators, respectively. The algebra su(l, 1) 
admits an infinite-dimensional irreducible representation 
where actions of K + , K_ and K 3 on basis 10 are 



K+ | n) = v/(n+ l)(n + 3) | n + 1} , 
K_\n) = ^/n(n + 2) \ n - 1) , 
K 3 \n) = (2n + 3) | n) . 

With basis on the algebra su(l , 1), the Hamiltonian Q 
is rewritten as H — (ti 2 /2ma 2 )K + K^, and a number 
operator N is given by 



N =(l/2)(K 3 -3). 



(3) 



Due to the dynamical algebra associated, the Barut- 
Girardello coherent states \z) , z G C, for the ISW are 
eigenstates of annihilation operator K _ 1121 



E 



n). 



On the other hand, a Diophantine equation is of the 
following form 



D(x u . 



,x k ) 



0, 



(4) 



where D is a polynomial with integer coefficients. In 
present terminology, Hilbert's tenth problem may be 
paraphrased as: Given a Diophantine equation of type 
0, we should build a procedure to determine whether or 
not this equation has a solution in non-negative integers. 
From the concluding results obtained by Matiyasevich, 
Davis, Robinson, and Putnam we know that, in the gen- 
eral case, this problem is algorithmically insolvable or 
more precisely, it is TM incomputable |6 1. 

With the basis for the ISW, our hypercomputational 
quantum algorithm for the Hilbert's tenth problem 1 5 1 
is shown in Table Q. Our algorithm is probabilistic 
such are quantum algorithms in general. Unlike these 
however, our algorithm is of adiabatic quantum compu- 
tation over infinite-dimensional spaces, given that part- 
ing from initial ground state (|5} of (|6), the final ground 
state J10I of 0, is obtained via the adiabatic theorem 
for unbounded operators 1 10 1. As we have stated else- 
where 1 5 1, the construction of our algorithm, breaking 
from Kieu's algorithm, opens the possibility to con- 
struct similar algorithms based on physical references 
whose associated dynamic algebra possesses an infinite- 
dimensional irreducible representation. 

3 Simulation Procedure 

Once constructed with (0 and (|6j selected for a 
Diophantine equation of type (0}, our problem revolves 
around solving (|8) in search of a value of T which sat- 
isfies condition (lilt , whereby we can establish criteria 
d!2i . Given that (|8j is formulated inside an infinite di- 
mensional Fock space, the simulation consists in being 
solved numerically in a Fock space truncated from di- 
mension m + 1, selected such that 1131 



m+l 



E 



has a norm less than one by some chosen e 



t \ 7 
\n) 



(z | z) 



m+l 



< e, 



(13) 



where | z) m+1 represents the infinite-dimensional state 
| z) truncate to the m + l dimension. 



Table 1. Hypercomputational quantum algorithm. 



Given a Diophantine equation with k unknowns of type 
J4j, we provide the following quantum algorithm to decide 
whether this equation has any non-negative integer solution 
or not: 

1 . Construct a physical process in which a system initially 
starts with a direct product of k coherent states \if)(0)), 
and in which the system is subject to a time-dependent 
Hamiltonian Hx(t) over the time interval [0,T], for 
some time T, with the initial Hamiltonian Hi and the 
final Hamiltonian Ho, given by 





= ®l*>, 
1=1 


(5) 


H A (t) 


= (l-t/T)H l + (t/T)H D , 




Hi 


k 

= E ( K +* - **) ( K -* - *) . 
j=i 


(6) 


H d 


= (D(Ni,...,N k )f. 


(7) 



2. Measure through the time-dependent Schrodinger equa- 
tion with the 

id t | V(*)> = H A (t) | V(*)> , for t € [0, T] (8) 

the maximum probability to find the system in a partic- 
ular number state at the chosen time T, 



where 



<(T) = max \(ip(T) \ m, 

(ni,...,n fc )SN fc 

= \ (i>(T) I ni,.. . ,n k ) \ 2 , 
\{n}) = |m, ...,n k ) - 



,nk)Y 



(9) 
(10) 



3. If -Pmax(T) < 1/2, increase T and repeat all the steps 
above. 

4. If 

Pmax(T) > 1/2 (11) 

then | {n}) is the ground state of Ho (assuming no 
degeneracy) and we can terminate the algorithm and de- 
duce a conclusion from the fact that 

Ho | {n}) = iff J4] has a solution in N. (12) 



The numerical solution involves a discretization pro- 
cedure that consists in replacing (|8} by a discrete equa- 
tion, whose solver preserves the norm of the state vector 
in the evolution of time. This solver is unitary and it is 
given by the Cayley transform of H™ +1 1131 



|V(*+1)> 



m+1 



1 - 



l rrriL 
2 H A 



1 



i Tjrn- 
2 H A 



my 



(14) 



4 Simulation parameters and results 

For the numerical simulations we have selected three 
Diophantine equations. Since we counted on all the 
information available for the simulation, the maximum 
probabilities (|9j can be calculated directly. Of course, 
this situation would be different if we were measuring 
upon the physical reference, in which case these prob- 
abilities should be obtained parting from the expected 
values of the operator numbers (|3j. 

4.1 Equation with solution 

We now consider a simple example which neverthe- 
less has all the interesting ingredients typical for a gen- 
eral simulation of our hypercomputing quantum algo- 
rithm. Let D(x) be the Diophantine equation 

D(x) = x - 6 = 0. (15) 

In relation to the values of z, these should be greater 
than 1.6 to guarantee that the maximum peak of density 
of probability associated with the initial ground state 
is less than 1/2 |5|. We choose z = 4 and m = 9 which 
satisfy Jl 3I > for e = 1CP 4 . The truncated Fock space has 
only ten dimensions, and is generated by { | 0) , . . . , | 9) }. 
ForO, 

according to the algorithm pointed out in Table 
0, we have that 

|^(0)} 10 « 

[0.16 0.36 0.51 0.53 0.43 0.29 0.17 0.08 0.04 0.02] T , 
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Applying dl4l > and incrementing T, figure Q in- 
dicates that the maximum probability reaches the 
value of 1/2 for T « 35 and then the Fock state 
| 6) 10 = [0 ... 1 0] T , is identified as the fi- 
nal ground state d 1 01 according to Jl li . Now, given that 



0, we can conclude that dl5> is solved 



#10X10 | g) 10 

in the non-negative integers in agreement with d!2i . 

4.2 Equation without solution 

Let D(x) be the Diophantine Equation 

D{x) = x + 6 = 0, (16) 




Figure 1. Simulation results for ITsl with (z, m) = (4, 9). 



and let z = 4 and m = 9. The figure (0 shows the prob- 
abilities of states as a function of T. Below T k 70 none 
of the states have probabilities greater than 1/2, and in 
fact the excited states, clearly dominate in this regime. 
Eventually, we enter the quantum adiabatic regime upon 
which the dominant state raises its probability over the 
1/2 value; indeed it corresponds to the Fock state | 0). 
Parting from dl 2i . it is concluded then that d!6i cannot 
be solved in the non-negative integers. 
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Figure 2. Simulation results for (TJJ with (z, m) = (4, 9). 



Unlike the previous example, this simulation does not 
give a conclusive result. In fact, it would be more pre- 
cise to affirm that according to figure (|2j, equation d!6i 
does not have a solution in {0, . . . , m = 9}. Therefore, 
it would be necessary to increase the value of m and ex- 
ecute the algorithm once again. Unfortunately, different 
to what would occur if we could count on the physical 
reference, in the case of the simulation, a standing cri- 
teria to determine the maximum value of m from which 
we can obtain a conclusive result, does not exist. 



4.3 Catalan Equation 

Given that every exponential Diophantine equation 
can be reduced to a polynomial Diophantine equation 
1 6 1, Hilbert's tenth problem also applies to these types of 
equations. It was recently proven that the Catalan Equa- 
tion x p — y q = 1, with x, y,p,q > 2, has one and only 
one solution (x = q = 3, y = p — 2) |14|. Although 
this solution could be found by simple inspection, differ- 
ent mathematical packages such as Mathematica™ or 
Maple™ , cannot find it. To put our algorithm to work in 
order to determine if the Catalan Equation can be solved 
in non-negative integers, we rewrite it as 

(a + 3) (b+2) - (c + 2)( <i+3 > = 1, (17) 

whose only solution is (a = b = c = d = 0). The 
simulation results obtained are shown in figure l|3}- Upon 
verification of d!2i . we can observe that our algorithm 
determines the existence of a solution. 
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Figure 3. Simulation results for {HJ wim ( z , m ) = (1-6, 3). 



5 Conclusions 

Although the numerical simulations cannot be equiv- 
alent to our possible hypercomputational quantum algo- 
rithm, due to the impossibility to simulate an infinite 
number of dimensions, their results give an idea of the 
behaviour that could be observed if we would have a 
real implementation and execution of our algorithm, sup- 
ported in a determined physical process. 

The nonequivalence is observed through an asym- 
metry between solvable and non-solvable Diophantine 
Equations in the non-negative integers, with respect to 
the conclusiveness of the numerical simulations. This 
asymmetry is usually a common characteristic of the 
Turing machine incomputable problems, for which, if 
the instance of the problem has a solution, a classic al- 
gorithm would find it; while if it doesn't have a solution, 
such algorithm will not halt. 
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